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We report on the direct observation of an oscillating atomic current in a one-dimensional 
array of Josephson junctions realized with an atomic Bose-Einstein condensate. The array 
is created by a laser standing- wave, with the condensates trapped in the valleys of the 
periodic potential and weakly coupled by the inter-well barriers. The coherence of multiple 
tunneling between adjacent wells is continuously probed by atomic interference. The square 
of the small-amplitude oscillation frequency is proportional to the microscopic tunneling rate 
of each condensate through the barriers, and provides a direct measurement of the Josephson 
critical current as a function of the intermediate barrier heights. Our superffuid array may 
allow investigation of phenomena so far inaccessible to superconducting Josephson junctions 
and lays a bridge between the condensate dynamics and the physics of discrete nonlinear 
media. 
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The existence of a Josephson current through a potential barrier between two supercon- 
ductors or between two superfluids is a direct manifestation of macroscopic quantum phase 
coherence (1,2). The first experimental evidence of a current-phase relation was observed in 
superconducting systems soon after the Josephson effect was proposed in 1962 (3), whereas 
verification in superfluid Helium has been presented only recently owing to the difficulty of 
creating weak links in a neutral quantum liquid (4,5). The experimental realization of Bose- 
Einstein condensates (BEC) of weakly interacting alkali atoms (6,7) has provided a route to 
study neutral superfluids in a controlled and tunable environment (8,9) and to implement 
novel geometries for the connection of several Josephson junctions so far unattainable in 
charged systems. The possibility of loading a BEC in a one-dimensional periodic poten- 
tial has allowed the observation of quantum phase effects on a macroscopic scale such as 
quantum interference (10) and the study of superfluidity on a local scale (11). 

A Josephson junction (JJ) is a simple device made of two coupled macroscopic quantum 
fluids (2). If the coupling is weak enough, an atomic mass current I flows across the two 
systems, driven by their relative phase A0 as: 



where I c is the "Josephson critical current", namely the maximal current allowed to flow 
through the junction. The relative phase dynamics, on the other hand, is sensitive to the 
external and internal forces driving the system: 



where h is Planck's constant divided by 2tt, t is time, and AV is the chemical potential 
difference between the two quantum fluids. Arrays of JJs are made of several simple junctions 
connected in various geometrical configurations. In the past decade, such systems have 
attracted much interest, because of their potential for studying quantum phase transitions 
in systems where the external parameters can be readily tuned (12). Recently, the creation 
of simple quantum-logic units and more complex quantum computer schemes (13) have 
been discussed. A great level of accuracy has been reached in the realization of two- and 
three-dimensional superconducting JJ arrays (12). One dimensional (ID) geometries are 
much more difficult to realize, due to the unavoidable presence of on site frustration charges 
that substantially modify the ideal phase diagram. ID JJ arrays with neutral superfluids 
(such as BEC), on the other hand, can be accurately tailored, and open the possibility to 
observe directly several remarkable phenomena not accessible to other systems (14)- First 
experiments with BECs held in a vertical optical lattice have shown the spatial and temporal 
coherence of condensate waves emitted at different heights of the gravitational field (10). 
More recently, the degree of phase coherence among different sites of the array ( 15) has been 
explored in the BEC ground state configuration. 

We report on the realization of a ID array of JJs by loading a BEC into an optical lattice 
potential generated by a standing wave laser field. The current-phase dynamics, driven by 
an external harmonic oscillator potential provided by an external magnetic field, maps on 
a pendulum-like equation and we have performed a measurement of the critical Josephson 
current as a function of the interwell potentials created by the light field. 

The experimental apparatus has been described in detail elsewhere (11). We produce 
BECs of 87 Rb atoms in the Zeeman state mp = — 1 of the hyperfine level F — 1 confined by 



I — I c sin A0 



(1) 



h— A0 = A^ 
dt 



(2) 
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a cylindrically symmetric harmonic magnetic trap and a blue detuned laser standing wave, 
superimposed on the axis of the magnetic trap. In essence, the cylindrical magnetic trap is 
divided into an array of disk shaped traps by the light standing wave. The axial and radial 
frequencies of the magnetic trap are, uo x = 2ir x 9 Hz and u r = 2tt x 92 Hz respectively. By 
varying the intensity of the superimposed laser beam (detuned 150 GHz to the blue of the Dl 
transition at A = 795 nm) up to 14 mW/mm 2 we can vary the interwell barrier energy Vo from 
to 5 Er where Er = h 2 /2m\ 2 is the recoil energy of an atom (of mass m) absorbing one of 
the lattice photons (16). The BEC is prepared by loading ~ 5 x 10 8 atoms in the magnetic 
trap and cooling the sample through radio-frequency-forced evaporation until a substantial 
fraction of condensed atoms is produced. We then switch on the laser standing-wave and 
continue the evaporation ramp until no thermal component is experimentally visible. This 
ensures that the system reaches the ground state of the combined trap. The BEC splits 
in the wells of the optical array: the distance between the wells is A/2 and ~ 200 wells 
are typically occupied, with ~ 1000 atoms in each well. The interwell barrier energy Vo, 
and therefore the tunneling rate, are controlled by varying the intensity of the laser, which 
is chosen to be much higher than the condensate chemical potential fi. fi ranges between 
fM 0.10 Vo for Vo = 2Er and fi ~ 0.04 Vo for Vo = 5Er. Each couple of condensates in 
neighbouring wells therefore realizes a bosonic JJ, with a critical current I c depending on 
the laser intensity. 

In a more formal way we can decompose the condensate order parameter that depends 
on position f and time t as a sum of wave functions localized in each well of the periodic 
potential (tight-binding approximation): 

*(r-t) = V^£^(*)$;(r) (3) 
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where Nt is the total number of atoms and ipj = yiijif) e*^ 3 '® is the j-th amplitude, with 
the fractional population rij = Nj/Nt and the number of particles Nj and the phase 4>j in 
the trap j. This assumption relies on the fact that the height of the interwell barriers is 
much higher than the chemical potential. We will prove by a variational calculation that 
this assumption is verified in most of the range of our experimental parameters (17). The 
wave function *&j(r) of the condensate in the j-th site of the array overlaps in the barrier 
region with the wavefunctions $j±i of the condensates in the neighbour sites. Therefore, the 
system realizes an array of weakly coupled condensates, whose equation of motion satisfies 
a discrete non-linear Schrodinger equation (18): 

iH l)r = ~ K ^ n ~ 1 + + ( £ ™ + A | 1p n \ 2 )lpn (4) 

where e n = Qn 2 , Q = |mu; 2 (|) 2 = 1.54 x 10~ 5 E R , and U = goN T J dr$j. The tunneling rate 

is proportional to K ~ — / dr [^;V$j • V$ J+ i + &jV ext &j + i}. A simple variational estimate, 
assuming a gaussian profile for the condensates in each trap, gives for Vo = 3 E R the values 
K ~ 0.07 Er, U ~ 12 Er and a chemical potential \x ~ 0.06 Vo that is much lower than the 
interwell potential Vo- We observe that the wavefunctions as well as K, depend on the 
height of the energy barrier Vo- 

Eq. (H) is a discrete non-linear Schrodinger equation (DNLSE) in a parabolic external 
potential, conserving both the Hamiltonian TC = Y,j{—K(4'j4'j+i + i'j'^j+i) + e j I i'j | 2 + j I 
ipj | 4 ] and the norm J2j n j — !■ 
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Although we can approximate the condensates in each lattice site as having their own 
wave functions, tunneling between adjacent wells locks all the different condensates in phase. 
As a result, when the condensates are released from the combined trap, they will show an 
interference pattern. This pattern consists of a central peak plus a symmetric comb of equally 
spaced peaks separated by ±2hkit exp / m where ki is the wave vector of the trapping laser and 
t exp is the expansion time. In practice one can think of the far field intensity distribution of 
a linear array of dipole antennas all emitting with the same phase. A complementary point 
of view is to regard the density distribution after expansion as the Fourier transform of the 
trapped one, i.e. the momentum distribution (19). It is easy to show that the sum of De 
Broglie waves corresponding to momentum states integer multiples of ±2hki is the sum of 
localized wavefunctions of Eq. (f|). The expanded cloud density distribution (Fig.|l|) consists 
of three distinct atomic clouds spaced by ~ 306 2hkit exp /m with the two external clouds 
corresponding to the first order interference peaks, each containing roughly 10% of the total 
number of atoms. The interference pattern therefore provides us with information about the 
relative phase of the different condensates (15,20); indeed, by repeating the experiment with 
thermal clouds, even with a temperature considerably lower than the interwell potential, we 
did not observe the interference pattern. 

This situation is different from the Bragg diffraction of a condensate released from a 
harmonic magnetic trap (21) where the condensate is diffracted by a laser standing- wave. 
In our case it is the ground state of the combined magnetic harmonic trap plus optical 
periodic potential that by expansion produces an interference pattern. For the time scales 
of our experiment the relative intensities of the three interference peaks do not depend on 
the time the atoms spend in the optical potential indicating that the steady state system 
has been reached. In absence of external perturbations the condensate remains in the state 
described by Eq. (^) with a lifetime of ~ 0.3 s at the maximum light power, limited by 
scattering of light from the laser standing- wave. 

In the ground state configuration the Bose-Einstein condensates are distributed among 
the sites at the bottom of the parabolic trap. If we suddenly displace the magnetic trap 
along the lattice axis by a small distance ~ 30 /im (the dimension of the array is ~ 100 /im) 
the cloud will be out of equilibrium and will start to move. As the potential energy that 
we give to the cloud is still smaller than the interwell barrier each condensate can move 
along the magnetic field only by tunneling through the barriers. A collective motion can 
only be established at the price of an overall phase coherence among the condensates. In 
other words, the relative phases among all adjacent sites should remain locked together in 
order to preserve the ordering of the collective motion. The locking of the relative phases 
will again show up in the expanded cloud interferogram. 

For displacements that are not very large, we observe a coherent collective oscillation of 
the condensates; i.e. we see the three peaks of the interferogram of the expanded condensates 
oscillating in phase thus showing that the quantum mechanical phase is maintained over the 
entire condensate (Fig||). In Fig. 2 A, we show the positions of the three peaks as a function 
of time spent in the combined trap after the displacement of the magnetic trap, compared 
with the motion of the condensate in the same displaced magnetic trap but in absence of the 
optical standing wave (we refer to this as "harmonic" oscillation). The motion performed 
by the center of mass of the condensate is an undamped oscillation at a substantially lower 
frequency than in the "harmonic" case. We will comment on this frequency shift later in the 
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text; we would like now to further stress the coherent nature of the oscillation. To do so we 
repeat the same experiment with a thermal cloud. In this case, although individual atoms 
are allowed to tunnel through the barriers, no macroscopic phase is present in the cloud and 
no motion of the center of mass should be observed. The center of mass positions of the 
thermal clouds are also reported in Fig.|2|B together with the "harmonic" oscillation of the 
same cloud in absence of the optical potential. As can be seen, the thermal cloud does not 
move from its original position in presence of the optical lattice. Indeed, if a mixed cloud is 
used only the condensate fraction starts to oscillate while the thermal component remains 
static, with the interaction of the two eventually leading to a damping of the condensate 
motion. 

We now turn back to the discussion of the frequency reduction observed in the oscillation 
of the pure condensate in presence of the optical lattice. The current flowing through the 
junction between two quantum fluids has a maximum value, the critical Josephson current 
I c , which is directly proportional to the tunneling rate. The existence of such a condition 
essentially limits the maximum velocity at which the condensate can flow through the in- 
terwell barriers and therefore reduces the frequency of the oscillations. As a consequence, 
we expect a dependence of the oscillation frequency on the optical potential through the 
tunneling rate. 

To formalize the above reasoning, we rewrite the DNLSE (Eq. |j) in terms of the canon- 
ically conjugate population/phase variables, therefore enlightening its equivalence with the 
Josephson equations for a one dimensional junction array: 

( hhj = 2K \JnfhjZ\ sin (<pj — 4>j-i) ~ 2K^/njrij + i sin {4>j+i ~~ 4>j) 

| h(j)j = —Urij — Vtj 2 + Kyjrij-i/nj cos (<f)j — + K^Jrij+i/rij cos (4>j+i — <j)j) 

It is useful to introduce collective coordinates (18): the center of mass and the 
dispersion a{t) are defined, respectively, as £(t) = J2jj n j an d cr 2 (t) = J2j j 2n j ~ £ 2 - From 
Eq. [5], we have h£ = 2K J2j y/ n j n j+i si n (<Aj+i ~ < Pj)- As the number of atoms is large, 
the "kinetic" energy term of DNLSE is small respect to the potential and nonlinear terms, 
and the population density profile is simply given by an inverted discrete parabolic profile, 
centered around £ (22): nj(t) = ^~ n C?~£) • furthermore J^cr 2 = 0. 

During the dynamical evolution, the relative phases across the junctions <pj+i — <pj = 
A0(t) remain locked together to the same (oscillating) value. This has been verified by 
numerically studying the Fourier transform ipk = Hj^fi %k ^\ from the experimental point of 
view this means that the expanded condensate continues to show the three peaks of the 
interferogram of Fig.Q. Therefore in these collective coordinate the current-phase relation is 
given by 

hf t at) = 2K sinA^t) . . 

hi**/,® = - mux* f m W 

which, in analogy with the case of a superconducting Josephson junction (in the resistively 
shunted junction model (1,2)) and with the case of 3 He (5), is a pendulum equation with 
the relative phase A(p corresponding to the angle to a vertical axis and the center of mass £ 
being the corresponding angular momentum. The current-phase dynamics does not depend 
explicitly on the interatomic interaction. This allows us to study regimes with a number 
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of condensate atoms spanning over different order of magnitude, which is different from 
the configuration considered in (10) where nonlinear effects would dephase the collective 
dynamics. However, it is clear that the nonlinear interaction is crucial to determining the 
superfluid nature of the coupled condensates, by locking the overall phase coherence against 
perturbations. 

From Eq. || we can see that the small amplitude oscillation frequency ui of the current 
/ = gives a direct measurement of the critical Josephson current I c = 2KNx/?i and, 

therefore, of the atomic tunneling rate of each condensate through the barriers. The critical 
current is related to the frequency uj of the atomic oscillations in the lattice and to the 
frequency u x of the condensate oscillations in absence of the periodic field by the relation 



mX 2 \uj x 

Figure ||] shows the experimental values of the oscillation frequencies together with the 
result of a variational calculation based on Eq. [5[ It must be noted that, due to mean 
field interactions, in our system a bound state exists in the lattice only for potentials higher 
than ~ Er] frequency shifts for lower potential heights are better explained in terms of the 
effective mass -^-jj = 0^ of the system (11). 

Increasing the initial angular momentum £o, the pendulum librations become anhar- 
monic, and can eventually reach the value A(f) max = tt/2. The system becomes dynamically 
unstable, and the phase coherence is lost after a transient time (the interference patterns 
washes out). In this regime, the pendulum analogy breaks down and a different dynamical 
picture would emerge. 

With this work we have verified that the BEC's dynamics on a lattice is governed by 
a discrete, non-linear, Schrodinger equation. This equation is common to a large class of 
discrete non-linear systems, including polarons, optical fibers, and biological molecules (23), 
thus opening up interdisciplinary research. The phase rigidity among different wells can be 
probed against thermal fluctuations to test various theories of decoherence ( 24 )■ One could 
study the role of collective dynamical modes in the creation of solitons and kinks of the 
type described in (18,25,26) (see also (23)) and more generally the routes to quantum phase 
transitions in nonhomogeneous, low- dimensional systems. 
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FIGURES 




FIG. 1. A) Combined potential of the optical lattice and the magnetic trap in the axial di- 
rection. The curvature of the magnetic potential is exaggerated by a factor of 100 for clarity. B) 
Absorption image of the BEC released from the combined trap. The expansion time was 26.5 ms 
and the optical potential height was 5Er. 
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FIG. 2. A) Center of mass positions of the three peaks in the interferogram of the expanded 
condensate as a function of the time spent in the combined trap after displacement of the magnetic 
field. Up and down triangles correspond to the first order peaks, filled circles correspond to the 
central peak. Open circles show the center of mass position of the BEC in absence of the optical 
lattice. The continuous lines are the fits to the data. B) Center of mass positions of the thermal 
cloud as a function of time spent in the displaced magnetic trap with the standing wave turned on 
(filled circles) and off (open circles). 
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FIG. 3. The frequency of the atomic current in the array of Josephson junctions as a function 
of the interwell potential height. Experimental data (circles) are compared to the calculated values 
(triangles). Each experimenatl data point was taken after a complete oscillation in the displaced 
magnetic trap. The oscillation was then fitted with a sine function giving the corresponding 
frequency (error bars are the standard deviation of the data from the fit. 
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